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An intcgrablc two-dimensional system related to certain fermion-soliton systems is studied. The self-consistent 
solutions of a static version of the system are obtained by using the tau function approach. The self-consistent 
solutions appear as an infinite number of topological sectors labeled by n € Z+, such that in each sector the 
scalar field would evolve continuously from a trivial configuration to the one with half integer topological charge. 
The spinor bound states are found analytically for each topological configuration of the background scalar field. 
The bound state energy satisfies an algebraic equation of degree 2n, so the study of the energy spectrum finds 
a connection to the realm of algebraic geometry. We provide explicit computations for the topological sectors 
n = 1, 2. Then, by monitoring the energy spectrum, including the energy flow of any level across E n = 0, we 
discuss the vacuum polarization induced by the soliton. It is shown that the equivalence between the Noether and 
topological currents and the fact that the coupling constant is related to the one of the Wess-Zumino-Novikov- 
Witten (WZNW) model imply the quantization of the spinor and topological charges. Moreover, we show that 
the soliton mass as a function of the boson mass agrees with the Skyrmes's phenomenological conjecture. Our 
analytical developments improve and generalize the recent numerical results in the literature performed for a 
closely related model by Shahkarami and Gousheh, JHEP06(2011)116. The construction of the bound states 
corresponding to the topological sectors n > 3 is briefly outlined. 
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I. INTRODUCTION 



Many aspects of soliton-fermion systems in two dimensions, such as the quantum corrections, have been 
studied in the literature; however, some important features at the zero order are yet to be investigated. In this 
context, the sl(2) conformal affine Toda model coupled to matter fields (CATM) is a two dimensional integrable 
and conformally invariant field theory with spinor and scalar fields, which presents topological solitons in the 
pseudoscalar field. It has two Dirac spinors coupled to a pseudoscalar field ip and a scalar field r\ which plays the 
role of a conformal symmetry 'gauge' connection. By choosing a particular constant solution r\ of the equations 
of motion the conformal symmetry is broken, thus obtaining a massive theory, the so-called affine Toda model 
coupled to matter (Dirac) fields (ATM) involving the Dirac spinors and the pseudoscalar field p only. The 
intcgrability of the theory is established using a zero curvature formulation of its equations of motion based on 
the si (2) affine Kac-Moody algebra. Some type of solutions such as the one and two soliton solutions, associated 
to some particular vacua, have been obtained in [1, 2]. Besides the conformal and local gauge symmetries, the 
model presents chiral symmetry and some discrete symmetries. However, one of its main properties is that for 
certain solutions there is an equivalence between the U(l) Noether current, involving the spinors only, and the 
topological current associated to the pseudoscalar soliton. That fact was established in [1] at the classical level 
and in [2] at the quantum mechanical level. It implies that the density of the U(l) charge has to be concentrated 
in the regions where the pseudoscalar field has non vanishing space derivative. The one-soliton of the theory is 
a kink type soliton with fractional topological charge, and therefore the charge density is concentrated inside 
the soliton. The CATM model constitutes an excellent laboratory to test ideas about confinement [2-4], the 
role of solitons in quantum field theories [2] and duality transformations interchanging solitons and particles 
[3, 5]. 

A version of the off-critical ATM model with just one Dirac spinor has been considered in many physical 
situations related to fermion-soliton systems. For example it appears in the study of the dynamical Peierls' 
energy gap generation in one-dimensional charge-density wave systems [6], quantum field-theory description of 
tunneling in the integer quantum Hall effect [7] , discussions of fractional charge or topological charges induced 
in the ground state of a Fermi system by its coupling to a chiral field [8-13] , and the spectrum and string tension 
of the low-energy effective Lagrangian of QCD 2 [4] . The model has also been considered in order to discuss the 
realization of chiral symmetries in (1 + 1) dimensions [14]. In the strong coupling limit it is closely related to 
the integrable variable mass sine-Gordon model presented in [15, 16]. Recently, a self-consistent solutions have 
been considered for a related version of the model with additional ip 4 potential using a numerical method [17] 
and the Casimir energy of that system has been computed, such that a static pseudoscalar field is prescribed 
with two adjustable parameters [18]. 

Recently, there has been significant growth in the study of fermion-soliton systems in the context of the 
localization and bound states of the fermion on the so-called braneworld models, following the original work 
of Rubakov and Shaposhnikov [19]. There are some models in which the brane Lagrangian considers the kink 
soliton of the ip 4 or sine-Gordon type systems in one dimension (sec e.g. [20]). Some related issues such as the 
particle transfer and fermion spectrum in braneworld collisions have also been considered [21, 22]. 
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In [2] the soliton type solutions of the ATM model associated to some type of vacua have been obtained 
using a hybrid of the dressing transformation and the Hirota's method. The basic idea was to look for vacuum 
configurations where the Lax potentials lie in an abelian (up to central terms) subalgebra of sl(2). The solitons 
are obtained by performing the dressing tranformations from those vacuum configurations. Such procedure led 
quite naturally to the definition of tau-functions, and then the Hirota's method was easily implemented too. 
The authors discussed the conditions for the solutions to be real, and evaluated the topological charges. The 
interactions of the solitons was studied by calculating their time delays. It has been shown that it is attractive, 
and the time delays have shown to be the same as those for the sine-Gordon solitons. An interesting aspect of 
the solutions is that when the two Dirac spinors of the theory are related by a reality condition, then either the 
soliton or the anti-soliton disappears from the spectrum. It has been interpreted as indicative of the existence 
of a duality involving the solitons and the spinor particles. In addition, the system can be mapped to either, the 
sine-Gordon model or the massive Thirring model through different gauge fixing procedures. In that scenario 
the solitons of the sine-Gordon model are interpreted as the spinor particles of the Thirring theory. Thus, the 
ATM theory is related to the sine-Gordon and the massive Thirring models in its strong/weak coupling sectors, 
respectively [3, 5]. 

An important property of the model under consideration is that it possesses special type of scalar and spinor 
fields in the intermediate field configurations in which they interact. In this scenario the role played by the local 
gauge symmetries on the physical soliton spectrum, as well as the Noether and topological currents equivalence 
deserve some clarifications. Here we show that the scalar field can represent topological solitons with infinitely 
many different shapes and the corresponding fermion becomes bound state solutions. The topological charge 
can take any real value as a consequence of the fact that the ATM model does not possess a particle conjugation 
symmetry [8] . In the intermediate region of field space in which the spinor and scalar fields interact to greater or 
lesser extent the system possesses rich physical properties. The choice of complicated background configurations 
associated to the scalar field renders the problem analytically unsolvable. Then one could resort to certain 
numerical solutions, in which case some important physical properties and concepts might not be clear. Some 
approaches have been proposed to circumvent this difficulty, they regard the scalar field as an external potential 
with prescribed shape, e.g. in [12, 13] it is assumed as an external infinitely thin (wide) soliton, whereas in [9] 
it is assumed to be a piecewise linear potential, allowing a set of analytic solutions for the spinor bound states 
[11]. 

In this paper we provide self-consistent solutions for a version of the ATM model with one spinor field 
coupled to the pseudoscalar, as considered in the literature. The static version of the system of equations 
allows a set of fermionic bound states corresponding to a family of soliton profiles for the scalar field. The self- 
consistent solutions appear as an infinite number of topological sectors labeled by an positive integer n G Z + , 
such that in each sector the scalar field would evolve continuously from a trivial configuration to the one with 
topological charge Q top = ±f • The bound states are found analytically for each topological configuration of the 
background field. The potential functional related to the scalar field is generated dynamically and its expansion 
would reproduce the scalar particle mass and $ 4 theory interaction terms. We obtain the energy levels of the 
bound states and investigate the relationship between the soliton, fermion and boson masses in each sector. 



4 



These relationships will be used to discuss the Skyrme conjecture regarding the ratio between the boson and 
soliton masses. Recently, in [17] the authors studied, by means of a numerical method, an analogous model by 
incorporating a scalar particle mass term and the $ 4 theory potential by hand. Our analytical calculations are 
compared to some of the numerical results of that reference, providing exact corrections to some quantities and 
identifying certain qualitative resemblances. 

We perform explicit computations of the bound state energies E n for the topological sectors labeled by 
n = 1, 2 as the background field evolves from a topological trivial configuration to a nontrivial one, for various 
soliton and spinor bound state parameters. The steps followed in order to construct the bound state energy E n 
corresponding to the sectors n > 3 will be outlined. Then, by monitoring the energy spectrum, including the 
energy flow of any level across E n = (n = 1, 2), we discuss the vacuum polarization induced by the soliton. 
The spinor charge quantization is discussed in view of the equivalence between the Noether and topological 
currents and the fact that the ATM coupling constant is related to the integer coupling constant of the WZNW 
model. 

The paper is organized as follows: section II summarizes the main properties of the ATM model. In section 
III a new set of fields arc defined: a scalar field $ associated to topological configurations and some spinor 
components £ a (a = 1,2,3,4) related to the fermion bound states. In section IV the self-consistent solutions 
are obtained through the tau function approach. In subsections IV A and IV B we perform explicit calculations 
corresponding to the topological sectors n = 1 and n = 2, respectively. The solitons and the corresponding 
spinor bound states are computed, as well as the bound state energies and the soliton masses in each case. The 
boson mass associated to the scalar $„ (n — 1,2) is computed providing the same formula for the both sectors 
n = 1,2. In the both cases the bound state energy E n is a solution of a relevant algebraic equation. In section V 
the vacuum polarization of the fermion induced by the soliton is discussed for any sector n £ Z + . In section VI 
we discuss some properties of the model, such as the currents equivalence and charge quantization. In section 
VII the ratio (soliton mass) /(boson mass) is discussed for the sectors n = 1,2. In the section VIII it is provided 
a brief discussion on the bound states related to the topological sectors n > 3 and in section IX it is presented 
the discussions of the main results. The paper presents 20 figures in order to illustrate the solutions and main 
results. The appendix A presents some results concerning the topological sector n = 2. 



II. THE MODEL 

Consider the two-dimensional field theory defined by the Lagrangian 

C = ~dwd»<p + i^d^ - ^ e W^75 ^ + Id^d^r] - I ml e 2 " (1) 

where <p, r\ and v are scalar fields, and ip is a Dirac spinor. It is defined r/5 = ip T 70, with tp being a second Dirac 
spinor. The sl(2) CATM Lagrangian (1) follows if one makes tp — > ^/| tp, v — > £ v, ip — > ip, tp — > ip in 
the model of [2]. The CATM model is related to the so-called two- loop Wcss-Zumino-Novikov-Witten (WZNW) 
model through a Hamiltonian reduction such that k = q £ Z, where q £ Z is the WZNW coupling constant. 
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In many applications of [2] it has been considered the reality condition 

4> = H ^ ( 2 ) 

where e$ is a real constant. It has been shown that the sign of plays an important role in determining 
the spectrum of soliton solutions. The theory (1) in [1] is defined for a general complex fields, i.e. a complex 
Lagrangian. In order to define a suitable physical Lagrangian one can follow the prescription to restrict the 
model to a subspace of well-behaved classical solutions [5] . For example, for = — 1 one can drop an overall 
minus sign from the Lagrangian in order to get usual kinetic terms for the spinor and <p scalar fields in such a 
way that one defines a submodel with a Hamiltonian bounded from below. In fact, as presented in [2] the one 
and two soliton solutions for ip and the corresponding spinor solutions satisfy the reality condition (2). 

One of the main conclusions in ref. [2] is that without the reality condition (2) the theory (1) has two 
Dirac spinors and it also has the soliton and anti-soliton solutions associated to tp, since both signs of the 
soliton charge are admissible. However, if one imposes the reality condition (2) the theory (1) looses one Dirac 
spinor and also one-soliton (two-soliton) solution, since for a given choice of only one topological charge 
Qtopoi = sign (e^p) \ [Qtopoi — sign (e^,) ] corresponding to one-soliton (two-soliton) or one-anti-soliton (two- 
anti-soliton), is permitted. That is indicative of the existence of some sort of duality between solitons and 
antisolitons associated to (p and the spinor particles ip, ip, respectively. So, one can have (without (2)) soliton 
and anti-soliton solutions with p> real. However, a charge zero solution corresponding to the scattering of a soliton 
and anti-soliton for p real docs not exist. Such solution exists, however, for p> complex but asymptotically real. 

The corresponding equations of motion are 

d 2 p = im^ J I ^ lb e^ +t ^l rr*1>, (3) 

V rZ 

i'fd^ = e-J-Vf if, (5) 

3 2 v = -^^v/*^- \m\e 2 \ (6) 

d 2 V = 0, • (7) 

The integrability properties, the construction of soliton solutions for i], p>, ip fields associated to a particular 
vacua, and many other properties were discussed in [I, 2]. We will be interested in special type of topological 
configurations associated to a scalar field $ and their relevant fermionic bound states \a- These new fields 
will be related to the original <p and ip,ip fields, respectively, by means of certain transformations. We start by 
reviewing the symmetries of (1) and the eqs. of motion (3)-(7). 

Left-right local symmetries. The Lagrangian (1) is invariant under the local U(1)l <£> U(1)r transformations 

/fc Fk 
-e R {x+) + \J-e L {x-) ■ v^v- (8) 

and 

e -» [( 1 +75)efl(a:+)-(l-75)ei,(a;-)] ^ . ^ _^ e « [(l+75)£H(x+)-(l-7 5 )ei,(x_)]^ ^ 



6 



Notice that the eqs. of motion are invariant under these transformations, whereas the Lagrangian (1) is 
invariant up to a surface term. 

U(l) global symmetry. Notice that, by taking sr (x + ) = —sl (ir_) = — \ 8, with 6 = const., one gets a global 
U(l) transformation 

p -> (p ; v -> v ; 77 ^ 77 ; tp ^ e l$ ip ; ip -> e~* e V (10) 

The corresponding Noether current is given by 

■7" =7^7"^, fl M J" = 0. 
fc 

Let us define the £7(1) charge as 

QtP = -j^ J dxJ°, 

where the normalization factor 1 /V has been introduced for later convenience. 

Chiral symmetry. In addition, if one takes sr (x + ) = El (#-) = — \ a, with a 
chiral symmetry 

V>^e 475 "V; ^ -> e -475 " if ; a ; z/ 1/ ; 77 ^ 77 (13) 

with the corresponding Noether current 

Jg = ^757^ + ; « = (14) 

Topological charge. One can shift the <p field as p — > + mir (ni G Z), keeping all the other fields 
unchanged, that the Lagrangian is left invariant. That means that the theory possesses an infinite number of 
vacua, and the topological charge 

Qtopoi. = J dxf , f = ^\[l^d„ p (15) 

can assume non trivial values. This model possesses fractional topological charges as we will see below. 
CP-like symmetry. In addition, the Lagrangian (1) is invariant under the transformation 

x + o x_ ; V ^70^ i ^ — ie'yoV' i V ^ V 3 > v ^ v (16) 

where e = ±1. The reality condition (2) breaks such CP symmetry, for any real value of the constant e^,. The 
CP symmetry is preserved if one takes = ±i. In addition, the theory (1) is not invariant under spatial parity. 
In general, solitons arc transformed into anti-solitons by parity transformations; so, it is very unlikely that the 
theory has the both type of solitons associated to <p for a given choice of in (2). However, maintaining the 
constraint (2) the both type of solitons arise for certain submodels defined for some gauge invariant fields, which 
turn out to be the sine-Gordon and the massive Thirring models, in the scalar and the spinor sectors of the 
ATM model, respectively [3, 5, 23, 24]. 



(11) 
(12) 

= const., one gets the global 
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Conformal symmetry. The model (1) is invariant under the conformal transformations (We are using x± = 
t ± x, and so, d± = \ (d t ± d x ), and <9 2 = <9 2 - d 2 x = Ad+d-.) 

x+ x+ = f(x+) , x-^x-=g(x-), (17) 

with / and g being analytic functions; and with the fields transforming as 

tp(x+ , x-) tp(x+ , x-) = if(x + , x-) , 

e -u{x+,x-) _^ e -f>(x+,x_) = yljt ^8 e - v { x+ , x _) ^ ^ 

e -»j(a; + _^ e -»Kx+ >*-) = (/' (g')^ e _77 ( x+ ' x ~^ , 

^(x+, af_) ^ </$(£+, x-) =e5( 1 +T 5 ) 1 °s(/')"'+5( 1 -75)iog(s')" i ^(a; + ,x_), 

where the conformal weight (5, associated to e _I/ , is arbitrary, and V transforms in the same way as tp. 
Associated to the conformal symmetry (17) there are two chiral currents 

J = -i\^ T \l+-y 5 )ip + yj^d + ip + d + ri, J = i^ T (l-^ 5 )i> + ^d_ip + d_r] (19) 

satisfying 

d-J = ; d+J = 0. (20) 

Observe, from (18), that the currents J and J have conformal weights (1,0) and (0, 1) respectively. Under 
the conformal transformations (17), the chiral currents transform as 

J{x + ) [Infix+T 1 (J(x + ) - [In /'(*+)]') , (21) 
J{x.) — ► [\ng\x.)]- 1 {J{x.) - [lng'(x-)}') (22) 

Then, given a solution of the model, one can always map it, under a conformal transformation, into a solution 
where 

J = ; J = 0. (23) 

Such a procedure amounts to gauging away the free field r\ [5]. In other words, (23) are constraints implement- 
ing a Hamiltonian reduction. So, for every regular solution 77 = 77 the CATM defined on a space-time (x + ,X-) 
will correspond to an off-critical submodel, the so-called affine Toda model coupled to matter (ATM). For the 
particular solution fj = the CATM and the ATM are defined on the same space-time. So, for 77 = n = constant 
the Lagrangian for the set of equations of motion corresponding to the Toda and Dirac fields becomes 

£ = -i<9^<^V + - e "° ^ etv/s/k v 75 V>- ( 24 ) 

This Lagrangian is not conformal invariant and defines the off-critical ATM model which we will consider in 
this work. One can easily check that the constraints (23), once we have set 77 — i] , are equivalent to 

1^1^=1^. (25) 
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Therefore, in the reduced model, the Noether current (11) is proportional to the topological current (15). 
That fact has profound consequences in the properties of such theory. For instance, it implies (taking np to be 
proportional to the complex conjugate of ip) that the charge density ip^ip is proportional to the space derivative 
of ip. Consequently, the Dirac field is confined to live in regions where the field p is not constant. The best 
example of that is the one-soliton solution of (1) which was calculated in [1] and it is given by 

ip = 2^^arctan (cxp (2m^ (x — x — vt) /y/l — (26) 

, 1/4 \ 

V 1 ""/ ,_ i( , 2 '»+^- i o-»')/Vi-» 2 / 



(27) 



1-ie m * 1 

v = - ^ log (l + exp (4m^, (x - x - vt) / \J\ - v 2 ^j ) - ^m^x + x_ (28) 
r, = (29) 

and the solution for ip is the complex conjugate of ip, implying = 1 in (2). Notice that, from (15), one has 
Qtopoi. = \ for the solution (26). In fact, this is a fractional charge, which is in accordance to the topological 
charge definition provided in [1]. 

Notice that the solution for ip is of the kink type soliton, and therefore d x ip is non- vanishing only in a region of 
size of the order of m^ 1 . In addition, the solution for ip resembles the massive Thirring model type soliton. One 
can check that (26)- (28) satisfy (25), and so is a solution of the reduced model. However, one must emphasize 
that these solutions are not truly sine-Gordon/massive Thirring solitons, e.g. its topological charge is +1/2 
whereas the 1— soliton of the usual sine-Gordon model possesses integer charge +1. 

We point out that the condition (25) together with the equations of motion for the Dirac spinors (4)-(5) imply 
the equation of motion for ip, namely (3). Therefore in the reduced model, defined by the constraints (23), one 
can replace a second order differential equation, i.e. (3), by two first order equations, i.e. (25). 

The model (24) with additional mass and ip 4 interaction terms has recently been considered in [17, 18], and 
it has been solved by means of a numerical method and the bound states of the fermion and the shape of 
the soliton have been obtained. In this type of model the pseudoscalar field boundary condition is a priori 
prescribed since the zeros of the potential provides the values the ip field gets at spatial infinity, in order to 
achieve finite energy and localized energy density solitons. Whereas, in the ATM model the scalar field soliton 
acquires, apart from riiir, new boundary conditions compatible with the dynamics, as will be seen below. The 
topologically nontrivial configuration of the scalar field <p of the ATM model, which is coupled to the ATM 
spinor field, can be thought to evolve continuously and slowly from a topologically trivial configuration to get 
a soliton configuration by a continuous gauge transformation, and the corresponding spinor bound states get 
a continuous set of configurations. Then the topological charge of the intermediate solitons associated to the 
field <p can get any real value. This is reminiscent of the ideas present in the so-called smooth bosonization 
method [25] in which the bosonic and fermionic formulations of the same theory is understood as two different 
gauge fixings of a larger gauge-invariant action containing the both fermionic and bosonic fields. Equivalence 
of the both sectors of the theory then amount to the usual gauge-fixing independence of the S-matrix elements 
in gauge theories. This method has been used to derive the bosonization rules for certain models. It implies 
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that the standard bosonization and fermionization prescriptions interpolate two extreme field representations 
of a "larger" theory which in a completely smooth way can be brought to almost any desired form in between. 
The equivalence between certain bosonic and fermionic theories (e.g. sine- Gordon/massive Thirring) then has 
profound consequences in the properties of the original gauge invariant theory, extending to a continuum of 
theories in which fermions and bosons appear interactively to a greater or lesser extent. 



III. TOPOLOGICAL CONFIGURATIONS COUPLED TO THE SPINORS & 

In the next steps we rewrite the eqs. of motion of the ATM model in terms of new real scalar fields $, Ci, C2 
and real spinor fields 6 (a = 1, 4) defined by the relationships 

^ = e-^/ 2 (a+*6), ^ = e iCl/2 (6+*6), (30) 

$ = y87fc^ + (Ci + c 2 )/2 (31) 

Our aim is to investigate certain intermediate configurations coupling the scalar <J> and the fermion field, such 
that any desired nontrivial topological configurations can be achieved by the field $. It is expected that the 
spinors 6 and the scalar $ will be the analogs of the massive Thirring and the sine-Gordon fields, respectively. 
So, the system (3)- (5), provided that 77 = 770 = const., becomes 



25+6 + d+( 2 6 - Mi 6 sin $ + Mi £3 cos $ = 0, (32) 

20+6 - d+( 2 6 + Mi £3 sin $ + Mi £4 cos $ = 0, (33) 

2<9_6 - d-d 6 - Mi £2 sin $ - Mi 6 cos $ = 0, (34) 

2<9_6 + d-Ci 6 + MiCi sin $ - Mi 6 cos $ = 0, (35) 
- [(66 + 6261) cos $ - (£i£ 4 - 66) sin $ - 



5 2 <f +^ 2 Ci + ^ 2 6] =^[*], Mi = W e*. (36) 

The appearance of the term involving £/i[<J>] in the eq. (36) deserves some clarifications. First, since the model 
possesses a local gauge invariance one needs to introduce a certain gauge fixing condition, so this term would 
play this role. Second, this term will be part of the potential for the scalar $, providing the vacuum points 
in order to find localized soliton solutions as in [17]. Third, in the static version of the eq. (36), the interplay 
between the functional i7i[$] and the combination (6 + 6) wm reproduce an equation invariant under parity, 
as will be seen below. 

The static version of the system (32)-(36) written in the form of an eigcnfunction equation for the Hamiltonian, 
with relevant mass and interaction terms for a $ 4 model, has been solved in [17] by numerical methods. In fact, 
setting 6 = E x_ + h + (x + ), (2 = —Ex + + h-(x-) and taking the static version of the system (32)-(36) one 
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gets 



£1 - E £ 2 - Mi £4 sin $ + Mi £3 cos $ 
£2 + £ 6 + Mi £3 sin $ + A#i£ 4 cos $ 
£ 3 + E 1 £4 + Mi £ 2 sin $ + Mi£i cos $ 
£4 — -E 1 £3 — Mi£i sin $ + Mi £ 2 cos $ 



0. 



0. 



0. 



(37) 



(38) 



(39) 



(40) 



I6M1 r 



(£16 + 6&) cos $ - (£i£ 4 - £26) sin $ + $" - (C) 



(41) 



fc L 



where E is the eigenvalue in the spinor sector and £ = (Ci +C2)/2- The role played by the auxiliary field £ will be 
discussed below in more detail. In the static version of the ATM system (37)-(41) the functional t/i[$] must be 
chosen such that, for every set of solutions {$, £ a , (a = 1, 2, ..4)} satisfying the system (37)-(40), the kink would 
remain as a self-consistent solution of the cq. (41). It is useful to consider the spinor components £ a [$] and 
the auxiliary field to depend functionally on <f> as discussed in [26] regarding the self-consistent solutions 
of the $ 4 theory coupled to a fermion. Recently, the problem of a self-consistent solution for a fermion coupled 
to static scalar field in the form of a kink (domain wall) [17, 27] has been considered mainly using numerical 
methods. The self-consistent approach to find fermionic bound states must be contrasted to the results in which 
the scalar field kink is considered as an external field as presented in [28] for the case of the ip 4 theory and in 
[29] for the case of the kink-antikink in the sine Gordon model. However, as we shall demonstrate in the present 
work, the presence of the fermion changes drastically the kink profile and the relevant spectrum. 

The mass and interaction terms of the $ 4 theory can be considered as the first terms in the expansion of the 
functional C/i[$] of our model. Therefore, the system (37)-(41) with = fc/16 and £ = 0, considering the mass 
and <I> 4 interaction terms in the expansion of E/i[$], will be equivalent to the model recently presented in [17]. 

In this paper we provide some analytical solutions to the systems of eqs. (37)-(41). In order to tackle this 
problem we follow closely the tau function approach developed for the ATM model [2]. This approach will be 
adapted conveniently for the static version of the ATM model (37)-(41) and implemented taking into account a 
variety of topologically nontrivial configurations for the scalar field $ and the relevant fermionic bound states 
associated to the field £ a . We are interested in solutions for the field $ representing a kink type solution such 
that 



We consider those solutions to the system (37)-(41) in such a way that the Hamiltonian is invariant under par- 
ity transformation, then the terms involving the auxiliary field ( and the functional U\ [<&] must take convenient 
forms in order to leave invariant the equation (41) under parity transformation. 

Under the above assumptions the solutions to the system (37)-(41) will be eigenfunctions of the parity 
operator. Then one must have 



$(x) 



(42) 



$(±00) = ±(p , $(0) = 0, 4> — const.. 



(43) 



£i(x) = -a£ 4 {-x), &{x) = a&i-x), o- = ±l, 



(44) 
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where a is the parity operator eigenvalue. Moreover, one imposes the following conditions to the spinors 

£„(±oo)->0, a=l,...,4; £1 (0) = £o = const. (45) 

An additional consistency condition must be satisfied by the solutions; in fact, from the equations (37)-(40) 
one gets the relationship 

& + £2 - £3 - & = ci = const. (46) 

Since all the spinor components vanish at x — > ±00, according to the boundary condition (45), from now on 
we will set c\ — 0. The above set of conditions will restrict considerably the possible forms of the auxiliary field 
C and the functional C/i [<&] as discussed below. In addition to the above set of conditions we require a finite 
energy and localized energy density associated to the field <J> such that it becomes a kink type solution of the 
system, so the functional Lq[$] in (41) must have a degenerate vacua. 



IV. SELF-CONSISTENT SOLUTIONS AND THE TAU FUNCTION APPROACH 

Next, we solve the system of equations (37)-(41) analytically and determine the eigenvalue E. In order to do 
this we will follow the tau function approach and assume that the expressions of the fields £ a and $ in terms of 
certain tau functions to be closely related to the ones which have been proposed in the relativistic ATM model 
[1, 2]. So, the pseudoscalar and the spinor fields of the model are assumed to depend on the tau functions 
{r[ n \ 4 n \ r£„, tU by the Ansatz 

T (") 

e »„/ n = e -ifl„^_ ) „= 1,2,3,... (47) 
= ["Li" + [-Ll" ^ ( 48 ) 



= -i{r|,„ [-^ ]" ^ 2 '" - 5f,„ [4) 1" ^'"i (49) 



= i{rt, n [-^\ n e- m * n - Tl, n [~^] n e ie ^}, (51) 



1 

where 8 n , Q a ,n are real parameters. Notice that we have incorporated a positive integer n € Z + in the expressions 
for the fields, each n labels a particular solution corresponding to a kink profile with specific nontrivial topological 
configuration $„ and a set of relevant spinor solutions £ 0j „. 

In order to find a self-consistent solution to the system of equations (37)-(41) we implement the procedure 
in two steps for each topological sector n = 1,2. First, we substitute the Ansatz (47)-(51) into the system 
(37)-(40) for the spinors coupled to the field $„ and determine the eigenvalue E n . Second, the eq. (41) is solved 
by requiring a self-consistency condition in each topological sector. In fact, in the static version of the ATM 
system (37)-(41) the functionals Ui[$] associated to each topological sector (n) must be chosen such that, 
for every set of solutions {$ n , £ aj „, (a = 1, 2, ..4)} satisfying the system (37)-(40), the kink $„ would remain as 
a self-consistent solution of the eq. (41). 
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We will concentrate below on some properties in the cases n = 1, 2 by constructing the explicit solutions and 
computing explicitly some of the related quantities, such as the topological charges, bound state energies, scalar 
field potentials and their vacua, kink and boson masses, as well as the kink and the corresponding fermionic 
bound state profiles, respectively. In section VIII we will provide a short discussion on the cases n > 3. 



A. The sector n = 1 and the topological charge Qtopoi = -f- 

Let us set n = 1 in (47)-(51) and consider the tau functions 

(52) 
(53) 
(54) 

where b a , b a are complex constants and k, Q\ are real parameters. For simplicity let us choose a ,i = 6\ (a = 
1,2,. ..4). 

The imposition of the parity conditions (44) to the £ a ,i fields defined in (48)-(51) for the tau functions 
(52)-(54) provide the parameters relationships 





= 1 + e" 




T (1) 


= 1 - 


- e i9l e 2KX 








= he KX , 




b 2 e KX , 


tL 


= b 3 e KX , 




= b A e KX 


tL 


= he KX , 


n. i - 






= he KX , 




= b A e KX 



&3 = —ia b 2 e 1 " , 63 = ia b 2 e l " , 64 = — icr 61 e 4 " , b 4 = iabie 
Consider the parametrization 

b 2 = R 2 e ta \ b 2 = P2 e-* a \ h = Pl e ia \ h = Pl e~ ia K 
The spinor components become 

L,l = (-i) a - 1 Pae^ x -^ +e ^ +(-D 



a-1 



,(2»0i) 



a = 1,2. 



(55) 



(56) 



(57) 



^ _|_ p (2kx— i6»i) v ' j _|_ e ( 2 nx+i9i) 

The remaining components £ 0)1 (a = 3,4) can be obtained from the above expressions (57) taking into account 
the parity condition relationships (44). Notice that the expressions (57) define real £ a .i functions. Taking into 
account the relations (48)-(51) and the condition (46) the parameters (56), must satisfy 



p\ sin (2ai — 6±) = p\ sin (2a 2 — 61). 



(58) 



Then, the simplest choice ot\ = a 2 , p 2 = ±Pi satisfies the above consistency condition (58). On the other 
hand, the conditions (42)-(43) applied to the function $1 as defined in (47) for n = 1 will restrict the possible 
values of the parameter Q\. So, the field $1 can be written as 

sin (0i ) e 2KX 



$1 (x) = -Oi + 2 arctan j — 



:}, 0ieI q = [-^ + 2nq, | + 2nq], 



(59) 



cos (61 ) e 2KX . 

where q e Z. The asymptotic values of the field are $i(±oo) = ±#1 (6\ G [— | , |]). The case n = 2 in (47) will 
provide a kink with asymptotic values $2 (±°°) € [— 7r , n] studied in the literature, see below. The angle 9\ is 
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-0.5/ 
-1.0 










^1.5 







FIG. 1. The graph of the kink corresponding to the parameters 8± = n/2, k, = 1.57. 

Notice that $ x (±oo) = ±n/2 



restricted by the conditions (42) and (43) to belong to the intervals I q defined for any q £ Z. The kink (59) is 
plotted in Fig. 1. Define the topological charge as 

y+oo 

itopol = 7^ 



3G 



= »n(+Q0)-$n(-Q0) 

2ir 

In the case at hand it becomes Q^ po i = ^, so it can take any non- integer value in the interval Q^ pol € 
[-1/2,1/2]. 

The region of validity of the angle 0\ must also be consistent with the relation (58) and the relationship 
between the parameters due to the normalization condition on the fermion sector. 
The boundary condition for the spinor £i is given by 

6 0=6 = Pi \ 2 9l ~ (62) 

The spinor is normalized as follows 

[hi\ 



+ °° r_L i 2 

= — x 



0! csc (6»i)(p? + pi) + p\ cos (2ai - 6 X ) - p\ cos (2a 2 - 6*i) (63) 
= 1. (64) 

In the Figs. 2 and 3 we present the graphs of the solutions £ Q , a = 1, ...4 defined in (48)-(51) considering 
the tau functions (52)-(54) for the both signs of the parity eigenvalues a = ±1, respectively. The parameters 
in these Figs, satisfy the condition (58) and the normalization condition (63)-(64). One can see that the kink 
parameters k, 6\ always depend on the spinor parameters through the eq. (58), the boundary condition (62) 
and the normalization condition (63)-(64), therefore the shape of the soliton always depends on the spinor 
parameters, and vice versa. 
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FIG. 2. The bound state of the fermion as a function of x for the positive parity a = 1. 
It is plotted for pi = p2 ~ 0.5, ai — 012 = —Si = —1.3, k = 1.35. Notice the spinor 
components relationships £1,1 (—as) = —£4,1(2;) and £2,1 {—x) = +£3,1(2;) 
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FIG. 3. The bound state of the fermion as a function of x for the negative parity a = — 1. 
It is plotted for p>2 = —pi = —0.5, ot\ = az = —Si = —1.3, k — 1.35. Notice the spinor 
components relationships £1,1 (—2;) = +£4,1(2;) and £2,1 (—2;) = —£3,1(2;) 



In order that the eq. (41) to be satisfied one must choose the functionals £ and U\ in a consistent way. The 
substitution of the relationships (47)-(51) given in terms of the tau functions (52)-(54) into the equation (41) 
will determine the forms of C and E/i[$i]. In fact, the eq. (41) is identically satisfied provided that 



c 







J7i[$i] = Aicos$i +A 2 cos(2$i) + ^ 3 cos(3$i) 
Actk 2 



A x = 3i?i 



4 cos 9i ( ■ 



Mo 



A, 



Bi = p\ cos (2a 2 -61) - p\ cos (2ai — 6\), B 2 



A 



4<JK 2 



B x +B 2 



A, 



M 2 

- p\ - 2pxp 2 cos (ai - a 2 ), 



R\B 2 



(65) 
(66) 



M 2 



Ri 



a M 7 



k ' 12(sin6>i) 2 ' 

Some comments are in order here. First, the functionals £ and U\ were chosen such that the eq. (41) is 
identically satisfied, and their properties in general depend on the scalar $1, and consequently their behaviors 
under parity transformation will depend on the form of $1. Taking into account the relationships (42) and (44) 
for the scalar and the spinors, respectively, under parity transformation, and in order to maintain invariant the 



15 



cq. (41) under the parity transformation one has conveniently identified the vanishing sector with (" = and 
the odd sector with the term involving the derivative t/([$i] = f^- In fact, f/i[$i] is an even function under 
$i — >• — $i as it is clear in (65). 

It is useful to consider the spinors £ [$i] as functional of the scalar field $i. So, the terms involving the 
spinors in eq. (41) will provide a functional — [/ 2 [<I>i]. Then one is able to write an equation for the scalar 

dU 

$ 1= — , 17= + U 2 [$i] (67) 

with 

k 2 

U\$i] = . . „ .„ [-4 cos 6>i cos $i + cos (2$i)l + const. (68) 

U ? = a^l \2 COS $ l{Pl + P\ + 6S 1 COS 6 l ~ 6 I/ 9 2 cos «2 COS (a 2 - 01 ) + 

b(sm #i) z 

p\ sinai sin (ai — 0i) — p\p 2 cos («i — a 2 ) cos0i] cos$ — 

4/91/92 cos (ai - a 2 )(cos$i) 2 + (p\ + p\) cos 2$i}. (69) 



Second, setting to zero the derivative of the potential 

C/'[$i] - (cos 0i -cos$i)sin$i = (70) 

one gets the vacuum points $i >va ,c = {01 + 2imi, 7rn 2 / rii, n 2 G Z}. The vacuum $i lt , ac = 0i was expected 
since it corresponds to the kink in (59). In addition, the vacua $i jt , ac = ^n 2 are inherited from the original 
ATM model in (24). A representative soliton interpolating these type of vacua is presented in (26). 
For the boson mass one gets 

2 _ d 2 i/[$i] , 
m = a$p kl=0 

= 2 K 2 (sec|) 2 . (71) 

The soliton mass associated to the field $i given in (59) is obtained from the relation M so i = J^(^i) 2 , so 
one has 

M sol = 4|k| cot (0i)]. (72) 

It is a remarkable fact that the energy density associated to the scalar in (67) [30] 

e(x) = \{&f + U[S> 1 ] = 2U[$ l ], (73) 

provides the soliton mass (72) for the kink (59). The soliton mass parameters k, Q\ will be related later to the 
spinor parameters through (58) and the spinor normalization condition (see below). 

The changes in the local features of the soliton, once the asymptotic values $i(±oo) are fixed, will be 
considered below to study the spectral flow, so for later convenience we compute the slope of the kink at the 
origin x = 

pi = 2k tan (0i /2) (74) 
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Next, we consider the problem of finding the bound state energy E\. Substituting the relationships (47)- 
(51) given in terms of the tau functions (52)-(54) into the equations (37)-(40) one gets a homogeneous linear 
system of equations for the unknown independent parameters {6i, bi\. Then imposing the condition of vanishing 
determinant for the 2x2 matrix formed by the coefficients in that linear system of equations, in order to find 
a non-trivial solution, one gets a quadratic algebraic equation for the bound state energy E\ 

8k 2 - Ml + Ml cos (20i )1 El + 2kM 2 sin (20i)^i + K 2 Ml - M? - K 2 Ml cos (20i) + 

Mf cos (20i ) - 2a Mi{[4k 2 - Ml + Ml cos (20i)]£i + kMI sin (20i)} = (75) 

where Mi = e^m^e^ . In Figs. 5, 6 we provide some contour-plots for _Eivs0i and in the Figs. 7 and 8 
we plot E\\sp\ for a special set of parameters, see below more on this point. Notice that E\{6\ = 0) = 
^(<T + eo)j e o = ±1 • So, -Ei (0i = 0) = a Mi, i.e., there is a positive-parity threshold bound state at E\ — M\ 
and a negative-parity threshold bound state at E\ = —Mi (see the bottom right graph in Fig. 6). 
In the infinite limit of the slope p, — > +oo (n — > +oo) one has 

£i(0i) = iaM 1 (l±|cos0i|) (76) 

whereas, in the limit \i — >• (k — >• 0) one gets 

E 1 (8 1 ) = aM 1 (77) 

The result (76) must be contrasted with the exact results in [11, 12] for the piecewise linear potential $ with 
a constant slope in the middle region. In this limit the authors have obtained E\ — M\ cos 9 (0 < 9 a < n). 
However, in those references the shape of the potential $ is prescribed. 

Let us investigate the boson mass dependence of the kink energy. Consider the expression (72) of the soliton 
mass and the boson mass (71), then one has 

M sol = {2V2[1 - 0x cot X ] cos (^)} to. (78) 

This is an exact result and shows that the soliton mass is directly proportional to the boson mass. However, 
its expansion in powers of 0i becomes M so i w -\/2to0 2 (1 — 9f/8 + ...). This can be compared to the <J> 4 theory 
exact result M*^ fe = ^ mQ\. 

The boson mass dependence of the bound state energy can be obtained substituting the relationship (71) for 
k into the eq. (75). In fact, this procedure provides an implicit eq. for E and to. In Fig. 10 it is plotted E vs 
to for Mi = 1. 

On the other hand, one has the ratio M so i/m = 2 for 0i = n/2. This ratio will be discussed below in 
connection to the Skyrme conjecture. In the Fig. 22 this ratio has been plotted in terms of 0i (see the curve 
labeled by n = 1). 



1. Special set of parameters 

Next we discuss a special set of parameters and the fcrmionic bound state energies in the case n = 1. For 
simplicity in what follows we set = a\, pi = p\. This choice satisfies identically the relationship (58), and 
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FIG. 4. The graph of the function aj ^- > for 9 X G 
[— 7r/2 , 7r/2] which must be equal to the ratio as 
in (80). 



the normalization condition (63) becomes 

r+oo r 4 , a 2 

j M [E^=^ CSC W ( 79 ) 

The fermion wave function normalization to unity provides a further restriction on the parameters, so 

^f = ^, (80) 
\k\ 0i 

In the Fig. 4 it is plotted the positive region of the function [ 81 g fll ](> 0) which appears in (80) defining the 
values of k once the parameter p\ is fixed. 

A general region of validity of 9\ must be consistent with the intervals Ig in (59), the relationship (58) and 
the normalization condition (63)-(64) determined for each choice of the set {a.2, ot\, p2, pi, k}. This amounts to 
generalize the simplest case (79)-(80) considered above and we will not pursue it further. 

Some comments are in order regarding the Figs. 5 and 6 corresponding to the energy spectrum of the n = 1 
sector. First, the graphs of the bound state energies for Q\ negative are the reflection of the positive part through 
the origin (0, 0) provided the parity labels are changed. Second, for Q\ positive, the number of bound states 
and their separation from the positive (negative) continuum increases as the parameter p\ is increased. Third, 
for some regions in the parameter space there are up to four bound states, two for each parity. 

In order to investigate the spectral flow dependence on the local features of the soliton, once the asymptotic 
values $(±oo) are fixed for a particular value of 6\, we have plotted E vs p± in Figs. 7 and 8. For example, 
one observes that the energy level crosses E = for 0\ = 7r/2 at p\ = 0.4 as shown in the top left of Fig. 7. 
Notice that in the limit 9\ — > one gets the Dirac sea corresponding to a trivial kink configuration as shown 
in the bottom right of Fig. 8. Moreover, from the Figs. 7 and 8 one notices that as Q\ decreases the positive 
parity bound states approach the positive continuum faster than the ones with negative parity do the negative 
continuum. 
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{E(0! ),<?!} [E(fli), 0i ] 

71 71 71 71 




FIG. 5. The fermionic bound state energies when M\ — 1 and the values pi — 
10,0.9,0.5,0.3. The signs ± indicate the parity of the bound states. Notice that for 
pi ( or k) — > +oo , the graphs look like the function (76) for the both parities. 



The Fig. 9 shows the total energy (soliton mass plus the fermionic bound sate energy with a = +1) as a 
function of the boson mass for Mi — 1 and some values of the angle 9\. Notice that for 0\ — it/ 2 one has a 
linear behavior resembling the result of ref. [17]. 



B. The sector n — 2 and the topological charge Qt opo i = -f- 

Here we closely follow the relevant steps of the previous construction. Set n = 2 and 8 a 2 = 2 (a = 1, 2, ...4) 
in (47)-(51) and consider the tau functions 

t (2) = 1 + e~ l ^e 2K \ t[ 2) = 1 + e^e 2KX , (81) 
i,2 = he™ + b n e 3KX , t£ 2 = b 2 e KX + 6 22 e 3KX , 

_3 t „kx i u „3k,x 4 l „kx i l „3kx foo\ 

r^ 2 = o 3 e + o 33 e , r^ 2 = o 4 e + o 44 e , (82) 

T ^2 = °l e + °ll e ) T £,2 = + 0226 , 

2 = o 3 e + o 33 e , r^ ?2 = o 4 e + o 44 e , (8dJ 
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FIG. 6. The fermionic bound state energies when Mi = 1 and the values pi — 
0.15,0.055,0.025,0.005. The signs ± indicate the parity of the bound states. Notice 
that for pi (or k) — > 0, they reproduces (77), for the both parities. 
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FIG. 7. The fermionic bound state energies when M\ — 1 and the values 9i — 
n/2, 7r/3, 7t/5, 7t/10. The signs ± indicate the parity of the bound states. Notice that 
for 9i — tt/2 and pi > 0.35 there are four bound states, two of them for each parity; 
whereas for pi < 0.35 there appears just one bound state for each parity. Notice that as 
8i decreases the number of positive parity bound states decreases. 
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FIG. 8. The fermionic bound state energies when Mi = 1 and the values 61 — 
7r/15, 7r/30, 7t/50, ty x 10~ 3 . The signs ± indicate the parity of the bound states. For 
0i — > the bottom right plot reproduces (77), for the both parities. 




FIG. 9. The total energy ei (soliton mass plus the fermionic bound sate energy with 
<r = +1) as a function of the boson mass for Mi = l,9i = vr/2, tt/3, 7t/4, n/10. Notice 
the linear behavior in the case n/2. So, it reproduces the result reported in [17] through 
a numerical approach. Notice that for a wide range of values of the parameters 8i and m 
one has ei > Mi and the linear behavior of ei. 



where b a , b a , b aa , b aai a — 1,2,3,4 are complex constants and k, 02 are real parameters. Then the parity 
condition (44) provides the parameters relationships 

h = -icrbte, hi = -iubi, b 2 = -iab 33 , b 2 2 = —i<rb$, b 3 = -icrb 2 2, (84) 

633 = —iob 2 , t>4 = —iffbu, &44 = — iab\. (85) 
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FIG. 10. The bound state energy of the fermion as a function of the boson mass m for 
Mi — 1 and 9\ — n/2. The ± signs indicate the corresponding parity. Notice that the 
curve E 1 (m) never approaches zero as m increases. 



Let us choose the next parametrization which is consistent with the relationships (84)-(85) 

b — n p lUa b — n e~ %CLa 



baa = 11ae ll3a , b aa = Tj a t 



-ifSa 



, a = 1,2, ...4. 



Pa = -da + rji = api, tj 2 = ap 3 , 773 = ap 2 , 774 = crpi- 

Then, taking into account the relations (48)-(51) and the condition (46) the parameters (86)- 
must satisfy 



p\ sin (2a?!) — p\ sin (2a 2 ) — p\ sin (203) + p\ sin (2a 4 ) — 2a P2P3 cos (a 2 — ct 3 







P1P4 cos (ai — 014) 

From (48)-(51) and (81)-(83), and considering the above parameterizations the spinors become 

.,„ ,, J g«(aa— 6 2 ) _|_ iuY a g[~*(7a+#2)] gxp (2kx) 



[1 +exp(2ra- i6 2 )] 2 
( 1 w-i e i (- tt °+ 9 ') - iaTg e ['(T.+92)] exp Qkx) \ 

\ L ) Mi , „-fl M2 J' a ~ i ' 2 



(86) 
(87) 
(88) 

for ci = 0, 



(89) 



(90) 



,+62)] exp 

[1 + exp(2Kx + i6 l 2 )] 2 

where Ti = pi/pi, T 2 = p 3 /p 2 ,^\ = «4, 72 = 0^3. The expressions (90) define real functions. The remaining 
components £ a ,2 (a = 3, 4) can be obtained from the above expressions (90) taking into account the parity 
condition relationships (44). 

From (47) for n = 2 and (81) the field $ 2 can be written as 



$2(3?) = -20 2 + 4arctan jy 



sin (9 2 ) e 2KX • 



hel q = [-l+2nq, I 



2irq] 



(91) 



cos (U 2 ) 

where q E Z. Notice that the asymptotic values of the field i»2(±co) = ±2#2 (#2 G [ — § , §])■ The angle 6*2 
is restricted by the conditions (42) and (43) to belong to the intervals I q defined for any q 6 Z. This kink is 
qualitatively similar to the kink however the kink in (91) approaches the asymptotic values (±202) G [ — T j 7r ]- 
These kinks are plotted in Fig. 11. So, according to the definition (60) the corresponding topological charge 
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■M-v) 



FIG. 11. The graphs of the kinks <3?„(:r) (n = 1,2) corresponding to the parameters Qi = 
6 2 = 7t/2, ki = 1.35, k 2 = 1.57 (p = 0.5). Notice that $i(±oo) = ±tt/2, $ 2 (±oo) = ±tt. 
At x — the slope of the kink 2 is greater than the one of the kink 1 for this choice of 
parameters. 

(2) 

becomes Q top = 26* 2 /7r, which can take any non-integer value in the interval g< — 1 , 1 > and the integer ±1 
charges just for 6> 2 — ±7r/2, respectively. The slope of the kink at the origin x = becomes 

/i 2 = 4k tan (6» 2 /2) (92) 

The region of validity I q of the angle 62 would be further restricted when the condition (89) on the parameters 
and the normalization condition on the fermion sector are imposed below. 
The soliton mass associated to the field <J> 2 given in (91) becomes 

M sol = 16|«| [l-0 2 cot (0 2 )]. (93) 

The potential associated to the scalar field $2 will be discussed below. 

1. Special set 0} parameters 

Following the same steps as in the case n = 1, one has that the eq. (41) is satisfied provided that the 
functionals £[^2] and U{ [$2] are chosen in a consistent way. The substitution of the relationships (47)-(51) for 
n = 2, given in terms of the tau functions (81)-(83) into the equation (41) will determine the expressions for £ and 
C/j [$]• For simplicity we perform this process for a special set of parameters p a = p Q , a a = a a , a = 1,2, ...4. 

— (2) ' dU^ 

So, one can identify the vanishing sector with £ = and the odd sector with U{ [$ 2 ] = g ^ for the functional 

(2) 

given by 

Ul 2) [$ 2 ] = 2(csc6» 2 ) 3 [Af 2 PoCos(2a ) + 4k 2 sin 6> 2 ] cos (^) [2cos6» 2 - cos(^)]. (94) 

(2) 

The expression (94) has been obtained by setting p a — p D , a a = a in the eq. (A2). In fact, U{ [$ 2 ] is an 

even function under $ 2 — > — $ 2 . Next, considering the spinors £ a [ < ! ) 2] as functionals of the scalar field $ 2 , the 

(2) ' 

terms involving the spinors in eq. (41) will provide {— XJ\ [$2]} associated to the functional 

U ( 2 2) [$2] = -2(csc e 2 f [M 2 pl cos (2a Q ) + 8k 2 sin 9 2 ] cos (-^) [2 cos 2 - cos (-^)] (95) 
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where the eq. (A3) has been used for p a = p a , a a = a D . Then one is able to write an equation for the scalar 

*Z = ^> UW=UW[* 2 ] + UP[* 2 ] (96) 

with 

4k 2 

[$1 = , , „ xo [-4 cos 9 2 cos ($2 /2) + cos $ 2 ] + consi . (97) 
(sm# 2 ) 

Setting to zero the derivative of the potential (97) 

U(2) ' m = J^^^ 0862 -^^ 8 ^^' (98) 

one gets the vacuum points $ va c — {26*2 + 27rni, 27m 2 / ni, n 2 € Z}. The vacuum $2,«ac = 202 was expected 
since it corresponds to the kink in (91). In addition, the vacua $2,uac = 27m 2 are inherited from the original 
ATM model in (24). For the boson mass one gets 

m = ~ on k=0 

= 2/c 2 (sccy) 2 . (99) 

Notice that the boson mass above has the same form as in the n = 1 sector (71). This was expected since 
these expressions are associated to the same elementary boson mass of the theory. On the other hand, the 
energy density associated to the scalar field in (96) e(x) = ^(Q'2) 2 + t^[ ( f>]^ 2 - 1 = 2[/( 2 )[$], provides the same 
soliton mass as in (93) for the kink (91). However, the soliton parameters will be related below to the spinor 
parameters according to the eq. (58), the boundary and normalization conditions. 

For later purpose let us consider the expression (93) of the soliton mass and the boson mass (99) and write 
the relationship 

M so i = {872[1 - 02 cot 2 ] cos ( — )} to. (100) 

This is an exact result and shows, as in (78) for the n = 1 sector, that the soliton mass is directly proportional 
to the boson mass. 

Next, we obtain the algebraic equation satisfied by the bound sate energy E 2 - Substituting the relationships 
(47)-(51) for n = 2, given in terms of the tau functions (81)-(83), in the equations (37)-(40) one gets a ho- 
mogeneous linear system of equations for the unknown independent parameters {61, b 2 , 63, 64}. Then imposing 
the condition of vanishing determinant for the 4x4 matrix formed by the coefficients in that linear system 
of equations one gets a quartic algebraic equation for the bound state energy E 2 , corresponding to the both 
parities u = ±1, respectively 



r±E\ + r 3 El + r 2 E\ + r x E 2 + r - uAkM^EI + Sl E 2 + s ) = 



(101) 
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where 



r 4 = -8Ml{sm6 2 ) 4 1 r 3 = 64kM^ cos 2 (sin 2 ) s 



r 2 = 128k 4 + 2M-f 



2M 2 - (12k 2 - 7 -Ml) cos (20 2 ) + (Hk 2 - 2M X 2 ) cos (40 2 ) 



l 



Mfcos (60 2 ) 



6k 2 + Ml + (6k 2 - 2M X 2 ) cos (20 2 ) + M\ cos (40 2 ) sin (20 2 ) 



r = 4Mf 7k 4 - 5k 2 M : 2 - M$ + (9k 4 + 2k 2 M\ + M 4 ) cos (20 2 ) + 3k 2 M 2 cos (40 2 ) (sin 2 ) 2 



s 2 



8k 2 + M\ - Ml cos (20 2 ) sin (20 2 ), si = 8k(cos 2 ) 2 2k 2 - M 2 + M 2 cos (20 2 ) 



s = Ml 



k 2 - M x 2 + (3k 2 + Ml) cos (20 2 ) sin (20 2 ). 



In the Figs. 17 and 18 we provide some contour-plots for E 2 vs6 2 and in the Fig. 19 we plot E 2 vsp a for the 
special set of parameters. 

The boundary condition for the spinor £1 is given by 



1, 



£i(0) = Co = ^[ sec (y)] 2 (Pi cos G-i + p 4 a sma 4 ) 



(102) 



The spinor normalization for a general set of parameters is provided in the appendix (Al). For a special 
choice of the parameters p a — p Ql a a — a D from (Al) one has 



J 

J — ( 



, ■'— J |K|(sm0 2 ) 3 



0=1 



|0 2 cr + sin (2a a ) sin0 2 — cos0 2 [0 2 sin (2a a ) + (rsin0 2 ]| 



= 1, 



In the particular case a a = from (103)-(104) one has 



4p 2 



(sin 2 ) E 



(103) 
(104) 

(105) 



\k\ 2 -isin(20 2 ) 

The Figs. 12 and 13 show the solutions £ aj2 , a = 1, ...4 defined in (48)-(51) (for n = 2) considering the 
tau functions (81)-(83) for a = ±1, respectively. The parameters in these Figs, satisfy the condition (89) and 
the normalization condition (103)-(104). We must emphasize that the kink parameters k, 2 always depend on 
the spinor parameters through the eqs. (89), the boundary condition (102) and the normalization condition 
(103)-(104), so the shape of the kink always depends on the spinor parameters and vice versa. 

In the Fig. 14 it is plotted the function g 2 J S i'" i ^ 2g2 - ) > which appears in (105) defining the regions of 
validity of 2 consistent with the normalization condition (103)-(104). 



V. VACUUM POLARIZATION 



The vacuum polarization of the fermion induced by the presence of the soliton has been considered in [8] 
by introducing a powerful method, called the adiabatic method. In this method, the topologically nontrivial 
configuration of the background scalar field which is coupled to the fermions is assumed to evolve continuously 
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FIG. 12. The bound state of the fermion £ a ,2 as a function of x for the positive parity 
a = 1. It is plotted for p — 0.5, a = 0, #2 = 0.32, k = 0.69 . Notice the spinor 
components relationships fi,2(— x) = —£4,2(2) and £2,2(— x) = +£3,2(2;) 




-6 -4 -2 2 4 6 -6 -4 -2 2 4 6 



FIG. 13. The bound state of the fermion £ a ,2 as a function of x for the negative parity 
a — — 1. It is plotted for p = 0.5, q = 0, 62 = 0.32, re = 0.69. Notice the spinor 
components relationships £i,2(— 2;) = +£4,2(20 and £2,2(— x) — —£3,2(2;) 



and slowly from a topologically trivial configuration. Later on in [12] the adiabaticity requirement was lifted such 
that the energy spectrum of the fermions (in first quantization) coupled to the background field is computed 
during the adiabatic process. Defining the fermionic charge of the soliton as charge of the ground state, 
and monitoring the energy spectrum, including the flow of any level across zero, they have deduced that the 
vacuum polarization induced by the soliton receives two contributions. First is the adiabatic contribution, 
predicted in [8], which is a linear decrease with slope —<p /7r, {<j) > 0) whenever the interaction is in the form 
^e 1 * 75 ^, [$(±00) = ±(f> ]. This is due to the occurrence of the spectral deficiency in the Dirac sea as 4> 
increases. The second is the non-adiabatic contribution which occurs only when one or more fermionic energy 
levels cross E = 0, since in that case the definition of ground state changes. In the case when one energy level 
crosses zero it must be considered a jump of +1 in the vacuum polarization since this energy level is filled after 
crossing E = in the vacuum state, by definition. The formula for the vacuum polarization (VP) for the type 



2G 

4p„ ! [sin(8 2 )] 3 




FIG. 14. The graph of the function in (105) for 62 G [—tt/2 , 7t/2] which must be equal to 
the ratio 

\k\ 

[E(m|. mj 

'5 10 
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-1 

5 10 

FIG. 15. The bound state energy of the fermion as a function of the boson mass m for 
Mi = 1 and 82 — > ir/2 for the both parities. The plot corresponds to either the positive or 
the negative parity. The curves never approach zero as m increases, so these curves have 
the same behavior as the ones for the case n = 1 (see Fig. 10). 




of interaction described above becomes [9, 11-13] 

$„(+oo) - $„(-oo) 

VP n = -i 



2tt 



1,2,. 



n8 v 
n 



(106) 



with possible jumps of +1 at configurations where level crossing occur. Notice that in the above eq. we have 
considered the VP n for each topological sector n = 1,2, of the model considered in this paper. 

In the Figs. 5 and 6 (case n = 1) and 17 and 18 (case n = 2) the energy spectrum of the fermions coupled 
to the background field is computed during the adiabatic process. For each plot it is possible to compute the 
vacuum polarization as the parameter 9 n increases using the formula (106) and considering the jump +1 every 
time a level crossing occurs. 



27 



JO 



15 



FIG. 16. The total energy e-z (soliton mass plus the fermionic bound sate energy) as a 
function of the boson mass m for Mi = 1, 0% — > 7r/2. Notice that the line departs slightly 
from the linear behavior observed in the n = 1 case (see Fig 9) in this limit. 

In the Figs. 7 and 8 (case n = 1) and 19 (n = 2) we have plotted the energy spectrum in terms of the spinor 
parameters p\ (n = 1), p Q (n = 2), respectively, for various fixed values <fr n (±oo) = n9 n (n = 1,2) the solitons 
may achieve asymptotically. These plots show qualitatively the role of the local features of the solitons, such as 
the change of the slope at the origin by varying the parameters p Q , p\, on the level crossing effects at E n = 0. 



VI. CURRENTS EQUIVALENCE AND CHARGE QUANTIZATION 

Finally, let us rewrite the currents equivalence (25) in terms of the scalar field $ and the spinors £ a ; so, taking 
into account (30)-(31) one can write the zero component of the eq. (25) as 

a— 1 

Taking into account £ = and the spinor normalization condition to unity for any n sector, the cases n = 1, 2 
are considered in (63)-(64) and (103)-(104), respectively, the integration of the relationship (107) provides 

n9 n = n = 1,2,3... (108) 

k 

In this way we have established the relation between the angle n and the parameters of the ATM model k, e^. 
Since k is related to the WZNW model coupling constant q £ Z by k = 4-, then from (108) one has the 
relationship 

^ = ^«g; n = l,2,... (109) 

Moreover, the VP n expression can get another formulation in terms of the parameter and the WZNW 
coupling q, i.e. VP q — — Therefore, it is a remarkable fact that the topological and Noether currents 
equivalence present in the ATM model implies that the parameter 9 n defining the asymptotic values of the 
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FIG. 17. The fermionic bound state energy versus theta angle (EVS82) as a contour 
plot for Mi = 1, p a = po, OL a = 0, (a = 1,2, ..4) and the values p — 5,1,0.6,0.5. The 
continuous and the dashed lines correspond to + and — parities, respectively. Notice that 
for p — > +00 there are some regions for 82 with four and two bound states, respectively, 
whereas for some regions there is none. 



scalar $(±00) is related, in the spinor sector, to the parameter which defines the reality condition imposed 
on the two spinors of the ATM model (2). 

The relationship (109) has some remarkable properties. First, since the parameters 8 n and define the 
topological and spinor charges through the eqs. (61), and (2) and (11)- (12) , respectively, together with (109) 
one can get 

Qtl = Q { :\ Qui = — > <2i 9) = — ■ (no) 
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FIG. 18. The fermionic bound state energy versus theta angle (E vs 62) as a contour 
plot for Mi — 1, p a — po, a a = 0, (a = 1,2, ..4), the values p — 0.4,0.3,0.2,0.1 and the 
regions for 9 2 = {[— w/2, tt/2], [-tt/3, 7r/3], [— 7r/4, 7r/4], [-tt/8, it/8]}, respectively . The 
continuous and the dashed lines correspond to + and — parities, respectively. Notice that 
as p decreases the regions for 62 with bound states become smaller and close to the origin, 
whereas the regions with zero bound states becomes larger. 



It is clear that Q top may assume positive or negative values depending on the sign of 6 n . Second, the relation- 
ship (109) shows that the ratio / 9 n must be quantized and that the U(l) charge for fixed value of q is quantized 
in units of the topological charge Q[l p , i.e. — [Q^J p ] n; Q$ p — ^; (0 O e [— tt/2, 7r/2]). Conversely, the 

topological charge Q^ p for fixed value of n may assume discrete values, i.e. Q^ p = [Q^] |; = 4e</,. 
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FIG. 19. The fermionic bound state energies when M 1 = 1, p a — po, a a — 0, (a — 1, 2, ..4) 
and for the values 2 = 0.41vr, 0.25tt, O.Itt, 0.033tt, 0.02tt, 1CT 3 x tt, respectively. The + and 
— parity signs are indicated close to the corresponding lines, respectively. Notice that as 
#2 decreases the positive (negative) parity bound states approaches the positive( negative) 
continuum and in general the number of bound states decreases. 



VII. RATIO (SOLITON MASS)/(BOSON MASS) 

Let us investigate the boson mass dependence of the kink energy in the cases n = 1,2. Then from the 
relationships (78) and (100) one has for n = ir/2, n = 1,2 the ratios M so i/m 

M so r I 2, n = 1 

1 (111) 

n = 2 

In the Fig. 22 this ratio has been plotted in terms of the angles n (n = 1,2). Notice that for the soliton 

with M S oi = 1 from (111) one has m = 0.125 in the n = 2 topological sector, which is ~ 8% smaller than 

the ratio between the experimental masses of the pion and the nucleon which is around pio n mass — _ q ^47_ 
1 r nucleon mass 

Therefore, if in the n = 2 sector of the ATM model the pseudoscalar field $(x) plays the role of the pions, 
solitonic excitations (with M so i = 1 for m = 0.125) could be considered to be nucleons. In fact, our exact result 
improves the numerical result in [17] which provided m numer , = 0.116, i.e. the numerical calculation result is 
about 7% less than the analytical result. 

For n = 2 case an implicit expression relating the bound state energy and the boson mass can be obtained 
substituting the relationship (99) for k into the eq. (101). In Fig. 15 the plot E vs m is shown for Mi = 1 in 
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FIG. 20. The ratio (kink mass)/(boson mass) in the interval 9\, 62 6 [— \ , £•]. 

the limit 82 — > tt/2 for the both parities. In that figure it can be seen that at m = 0.125 there is no fermionic 
level crossing zero. Therefore, in this case the soliton of mass M so i = 1 polarizes the vacuum and its ground 
state possesses fermion number one, the vacuum polarization being (—1) for 82 = 7r/2 according to the formula 
(106). 

It could be interesting to consider the extension of the ATM model [23, 24, 32-34] incorporating more pions 
and nucleons to mimic the Skyrme model in order to partly account for the difference with the experimental 
value of the above result regarding the ratio between the pion and nucleon. 



VIII. THE TOPOLOGICAL SECTORS n > 3 

Following the previous constructions one can extend to the case n > 3. The tau functions in (47)-(51) can be 
written as 

r (,l) = l + e- ie "e 2KX , r[ n) = [4 n) }* (112) 
rl n ^b a e KX + b aa e 3KX + ... + baa...ae^ n -^ KX \ a = 1,2,. .4 (113) 

"ti" times 

^a T. „ KX 1 T ?,kx 1 1 T \(2n— 1)kx] (MA\ 

r e „ = o Q e + b aa e + ... + b aa^a e u ' (114) 

Let us discuss the appearance of the algebraic equation for the eigenvalue E in the case n > 3, the analog 
of the quadratic (75) and quartic (101) equations in the cases n — 1,2, respectively. Since the eqs. (39)-(40) 
are related to the eqs. (37)-(38) by parity transformation it is enough to consider the last couple of eqs. in 
order to get the mentioned algebraic equation for the eigenvalues E. Substituting the relationships (47)-(51) 
given in terms of the tau functions (112)-(114) into the equations (37)-(38) will provide a 2n— order algebraic 
equation for the bound state energy E. In fact, notice that the parameter E appears as a coefficient of the 
fields £1 and £2 in (37)- (38), and this system of eqs. will provide a homogeneous linear system of eqs. for the 
unknown parameters b a , b aa , ...baaa^a, K, baa, ■■■baaa^a, (a = 1,2). Moreover, the parameters of type b aa ,,, a 

depend on the other set of parameters b aa a through the parity relationships (44) [in the cases n = 1,2 one 
has the eqs. (55) and (84)-(85), respectively ]. So, in view of the expressions of the eqs. (37)-(38) in terms of 
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the relevant tau functions, the eigenvalue parameter E will appear only in the coefficients of the independent 



parameters 



61, &n, ... 6 H...1 &2, &22, ---b 22... 2 • Then, in order to get a non-trivial solution one must impose 



a condition det[^4 p9 ] = 0, where the matrix components A pq (p, q = 1, 2n) are the coefficients of the 2n— order 
homogeneous linear system of eqs. for the above independent parameters. The outcome will be an algebraic 
equation of order 2n for E with coefficients depending on the parameter constants (k, Mi, a, n ,). 

Finally, if one considers the eqs. (75) and (101), respectively as plane curves of fourth and sixth degrees in 
the variables {E, k} we can imagine that the possible combinations of coefficients (depending on the parameters 
Mi,6 n ,a) in the relevant algebraic equations will give rise to various families of plane curves; in this way the 
problem of the spectral flow of the fermion-soliton model due to the variations of the local features of the soliton, 
considered here for each topological sector "n", will be connected to the study of the plane curves of 2(n + 1) 
degree in the realm of algebraic geometry (see e.g. [35]). Notice that k is related to the soliton and fermion 
bound state parameters through the normalization conditions like (63)-(64) and (103)-(104), and this parameter 
determines the local properties of the soliton such as the slope at the origin for each value of 6 n , as in (74) and 
(92). 



IX. DISCUSSION 

We have considered special type of fields in the affine Toda model coupled to matter providing the analytical 
solutions as bound states of the fermion and solitons of the scalar field. The solutions of a static version of the 
system are obtained by using the tau function approach which has been previously proposed in the relativistic 
version of the theory. The self-consistent solutions appear as an infinite number of topological sectors labeled by 
n e Z + , such that in each sector the scalar field would evolve continuously from a trivial configuration to the one 
with half integer topological charge ±|. The spinor bound states are found analytically for each topological 
configuration of the background scalar field. In these developments the soliton shape depends crucially on 
the spinor bound state parameters, i.e. it depends on the fermionic state to which it is coupled. It is worth 
mentioning that we have considered the back-reaction of the spinor on the soliton exactly by using our analytical 
methods, this may be contrasted to the previous results in the literature in which the spinor back-reaction has 
been neglected or either considered perturbatively. 

The bound state energy in each topological sector satisfies an algebraic equation of degree 2n with coefficients 
depending on the parameters k, Mi, a, 8 n , so the study of the energy spectrum finds a connection to the realm of 
algebraic geometry. We provided explicit computations for the topological sectors n = 1,2. Then, by monitoring 
the energy spectrum, including the energy flow of any level across E n = 0, we have provided the formula in 
order to compute the vacuum polarization induced by the soliton. In this connection we have displayed some 
figures in the both topological sectors and elaborated on the adiabatic and non-adiabatic contributions to the 
vacuum polarization. 

It is shown that the equivalence between the Noether and topological currents present in the ATM model 
and the fact that the coupling constant is related to the integer coupling of the Wess-Zumino-Novikov-Witten 
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(WZNW) model imply the spinor charge quantization. 

Considering the spinor fields as functionals of the scalar field we have found their contribution to the potential 
associated to the scalar field, thus computing a soliton mass and the relevant elementary boson mass. In the n = 
2 topological sector we have determined that the ratio m/M sol approaches (boson mass) /(soliton mass) w 0.125, 
an error of about 8% with respect to the 4D experimental result which is around (pion mass) / (nucleon mass) w 
0.147, in this way reproducing the conjecture of Skyrme. So, we have extended the results of ref. [17] con- 
cerning several properties of the bound states and improved the result regarding the Skyrme's conjecture. The 
construction of the bound states corresponding to the topological sectors n > 3 was briefly outlined. 

Finally, in [1] the authors discussed a linear dependence on the fermion charge of the both soliton and 
spinor particle masses, M soliton <~ Q^, <~ however, the solitons discussed there were of the type (26) 
associated to the vacua (p vac . — nni, ri\ £ Z, ip vac . — 0, and the spinor particle mass was defined by expanding 
the interaction term in (24) around that soliton vacua. It could be interesting to discuss such relationships 
in the framework of our analytical solutions, in a future investigation. We may argue that those type of 
relationships are concerned with the correspondence (spinor) particle/soliton in the ATM theory, indicating 
some sort of duality similar to the electromagnetic duality of some four-dimensional gauge theories possessing 
the Bogomolny (monopole) limit [31] , where the masses of particles and monopoles (dyons) are given by mass ~ 
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We present some results related to the topological sector n = 2. The normalization of the spinor function in 
terms of the full set of parameters becomes 




Appendix A: Spinor normalization and the functionals J/i^^n] for n = 2. 




(Al) 
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Next, we provide the general form of the functional [$2] and U 2 2 ^ [$2] in terms of the full set of parameters. 
Considering i?[$ 2 ] = gfgggj one has 

r e( 2i9 ')M 2 (e ifla + R[^2]) 2 [h -b 2 + (hi - b 22 )R[^ 2 ]][b s ~b 4 + (fe 33 - M^N 
I (l + e^i?[$ 2 ]) 6 

e( 2 ^)M 2 (l + e^i?[$ 2 ]) 2 [&i - b 2 + {bu - b 22 )R[^ 2 ]][b 3 - fe 4 + (633 - fe 44 )i?[$ 2 ]] 



+ 



+ 



+ 



and 



[/ 2 (2) [<J>2] 



(A2) 



e( 2 ^)M 2 [bi - 6 2 + (611 - £22)^2]] [fe 3 + &4 + (fe 33 + b M )R[® 2 }} 
(l + e^i?[$ 2 ]) 4 

e( 2 ^)M 2 [b! + fe 2 + (fell + &2 2 )^[$ 2 ]][&3 ~ h + (fe 33 - Wgtgg]] 

(l + e^i?[$ 2 ]) 4 

e(- 2 ^)M 2 [6! + b 2 + (fe n + fe 22 )i?[$ 2 ]][fe 3 - &4 + (fe 33 - Mglggj] 

(l + e(-^)i?[$ 2 ]) 4 
e (-2»e 2 ) M2 [ bi _ &2 + (bn _ fe 22 )j?[$ 2 ]][fe 3 + b 4 + (fe 33 + fe 4 4)fl[$2]] 

(l + e(-^)i?[$ 2 ]) 4 

M 2 [fei + b 2 + (fen + b22)J?[^ 2 ]][fe3 + b A + (fe 33 + fe44)A[$2]] 

|l + e( <fl »)i2[$ 2 ]| 4 " " 

M 2 [fei + fe 2 + (fell + fe 22 )J?[^2]][fe 3 + fe 4 + (fe 33 + fe^flfoj]] 

|l + eM»)i2[$ 2 ]| 4 
16ze(-^) [-1 + e( 2ie2 )]K 2 (csc [i(20 2 - $ 2 )]) 2 sin 2 sin ^) 
|l + e(^)#[$ 2 ]|4 , 

1 e (- 4l92 )M 2 i?[$2] 
: 2 |l + e(^)i?[$ 2 ]| 6 X 

{-(e< i0a > + i?[$ 2 ]) 8 [(fen - fe 22 + (fei - fe 2 )i?[*2])(fe 33 - fe 44 + (63 - b 4 )R[$2])} + 

[(fel - fe 2 + (fell - fe 2 2)i?[$2])(fe 3 - fe 4 + (fe 33 - b 44 )R{$ 2 })}(l + B^Rfa]) 8 - 

[(fei + fe 2 )(fe 3 + 64) - (fen + fe 22 )(fe 33 + fe44)]e (4l92) |l + e (l92) i?[$ 2 ]| 4 (-l + i?[<5 2 ] 2 ) - 
iae^\l + eS-^R{<S> 2 }\ 2 (l + e^i?[$ 2 ]) 4 [(feii + fe 22 + (fei + 62)^*2]) X 
(61 - fe 2 + (611 - b 22 )R{<P 2 }) + (633 + 644 + (6 3 + 6 4 )i?[*2])(fe 3 - fe 4 + (fe 33 - fe 44 ) x 
R[$ 2 })} + iae^ll + J2[* 2 ]| 2 (l + e^ 6 ^ R[$ 2 ]) 4 [{bu - fe 22 + (fei - 62) x 
i?[$ 2 ])(fei + fe 2 + (fen + b 22 )R[<5> 2 ]) + (fe 33 - fe 44 + (63 - b 4 )R[$ 2 })(b 3 + fe 4 + 

(fe 33 +fe 4 4)i?[$2])]} (A3) 
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